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ABSTRACT 


As a plausible explosion model for a Type I supernova, the evolution 
of carbon-oxygen white dwarfs accreting helium in binary systems has been 
investigated from the onset of accretion up to the point at which a 
thermonuclear explosion occurs. Although the accreted material has been 
assumed to be helium, our results should also be applicable to the more 
general case of accretion of hydrogen-rich material, since hydrogen shell 
burning leads to the development of a helium zone. 

Several cases with different accretion rates of helium and different 
initial masses of the white dwarf have been studied. The relationship 
between the conditions in the binary system and the triggering mechanism 
for the supernova explosion is discussed, especially for the cases with 
relatively slow accretion rate. 

It is found that the growth of a helium zone on the carbon-oxygen 
core leads to a supernova explosion which is triggered either by the off- 
center helium detonation for slow and intermediate accretion rates or 
by the carbon deflagration for slow and rapid accretior. rates. Both 
helium detonation and carbon deflagration are possible for the case of 
slow accretion, since in this case the initial mass of the white dwarf is 
an important peurameter for determining the mode of Ignition. Finally, 
various modes of building up the helium zone on the white dwarf, namely, 
direct transfer of helium from the ccxnpanion star and the various types 
and strength of the hydrogen shell flashes are discussed in some detail. 


Subj<ot h«adinq>> stars; accretion - stars: binaries - 

■tarst evolution - otars: interiors - stars: supernovae 
- stars: white dwarfs 




I . INTRODUCTIOH 


a) Type I Supeznovae 

The origin of Type I supernovae (SN I) has long been obscure (see 
Wheeler 1981 for a review) . The spectra of SN I show no strong hydrogen 
lines, which implies that the progenitor stars of SN I should be hydrogen- 
deficient stars (e.g. Oke and Searle 1974) .. SN I have been observed in 
elliptical galaxies as well as in spiral and irregular galaxies (e.g. 
Tammann 1974) . SN 1 are not concentrated in spiral arms (Maza and van 
den Bergh 1976) . These facts have led to the ideas that the progenitor 
stars of SN I are white dwarfs (Finzi and Wolf 1967; Whelan and Iben 1973) 
or helium stars (Wheeler 1978; Arnett 1979) . 

Another characteristic of SN I is the exponential tail in their 
light curves (e.g. Barbon, Ciatti, and Rosino 1973). Such a tail requires 
an additional energy source at late times. Several energy sources have 
been proposed (see Chevalier 1981a for a review) . Among them, the radio- 
active decay model (®*Ni ■+ *®Co -► ®®Pe) has recently made remarkable 
advances (Arnett 1979; Colgate, Petscheck, and Kriese 1980; Axelrod 1980 
a,b; Chevalier 1981b; Weaver, Axelrod, and Woosley 1980). It has been 
shown that energy deposition of y-tekya and positronr emitted from the 
decays of ®®Ni and ®®Co can produce a good fit to the observed light 
curves of SN I. This model is supported by the good agreement between 
the observed spectra of SN I 1972e at late times (Kirshner and oke 1975) 
and the calculated synthetic spectra of Fe (Meyerott 1980) ^uld of Fe with 
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Co (Axolrod 1980a, b). 

Tho above success gives an important constraint on the progenitor 
sodels for SN I, i.e., relatively large aroo'ont of ’*^Ni should be produced 
through the primary explosion. Therefore we need to investigate whether 
the candidate stars of SN l actually evolve to produce sufficient ’*'Ni 
during the explosions. I have made such a computation for the evolution 
of accreting white dwarfs from the beginning of the mass accretion through 
the explosion. 


b) Evolution of Accreting White Dwarfs 

A white dwarf in a close binary system can evolve into a supernova 
explosion in the following way. When the companion star begins to expand 
over its Roche lobe, the white dwarf accretes the matter transferred from 
the companion star. When a sufficient amount of hydrogen-rich material 
is accumulated on the white dwarf, hydrogen shell burning is ignited. 

The hydrogen burniiig can proceed either in a thermally stable fashion or 
through unstable shell-flashes. 

The strength of the shell-flashes depends mainly on the accretion 
rate (see Sugimoto and Miyaji 1980 for a review) , Rapid accretion leads 
to a stable hydrogen burning or a relatively weak shell flash which pro- 
cess most of the accreted matter into helium. Slower accretion results 
in a stronger flash (sometimes nova-like explosion) which ejects a part 
of the accreted matter into space while processing the rest of the 
material into helium. Such hydrogen shell flashes recur many times 
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during th« accretion, and a helium aone gradually grotvj in mase near the 
surface of the white dwarf; this is equivalent to the accretion of helium. 
It is also possible that helium is directly transferred from the conq>anion 
star because, in sane cases of mass exchange, the coiqpanion star might 
have evolved into a helium star before filling its R<x:he lobe. (See §v 
for a more detailed discuss' on on the relationship between the strength 
of the hydrogen shell flashes, the growth of the helium zone and the 
conditions in the binary system and for a discussion of the possible 
evolutionary scenario leading to the direct accretion of helium.) 

The effect of the growing helium zone on the evolution of accreting 
helium white dwarfs was first investigated by Nomoto and Sugimoto (1977) . 
They simulated the average evolution by conputing the steady accretion 
of helium, which is a good approximation as far as the phenomena in the 
relatively oeep layers interior to the hydrogen-burning shell are con- 
cerned. It was found that relatively rapid accretion ignites the wea)( 
off-center helium flash while slow accretion leads to the ignition of 
helium at the center. In the latter case, the white dwarf disrupts com- 
pletely as a helium detonation supernova and ejects a large amount of 

(see also Mazurek 1973). 

Recently the same approach was adopted by Pujimoto and Sugimoto 
(1979, 1981, hereafter referred to as PS) and Taam (1980a, b, hereafter 
referred to as TM) for the study of long term evolution of accreting 
carbon and oxygen (C-fO) white dwarfs. They showed that, in some cases 
of accretion, the building up of a subst^mtial helium zone on the C-fO 
white dwarf leads to a very explosive helium ihell flash that could 
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ttiggar m ■upernova-lik* «xpIoBion. However, their reeulte were preli- 
minary becauae they asaumefl hydroetatic equilibrium in their CQnq;>utations 
and stopped computing at the beginning of the hydrodynamical stages. 

c) Present Study 

I have investigated the evolution of C-fO white dwarfs by adopting 
the same method as Nomoto and Sugimoto (1977) slatting from the onset of 
accretion of helium through the thermonuclear explosions, i.e., from the 
hydrostatic through hydrodynamic stages. Several cases with different 
accretion rates of helium and different initial masses of th? white dwarf 
have been ccmiputed. 1 have found that the accreting C+0 white dwarfs 
actually evolve into supernovae of various types, neunely, double detona- 
tion, single detonation, and carbon deflagration, depending on the 
conditions at the initiation of the explosion. Most of these models 
produce a large amount of ®^Ni as required by the radioactive decay model 
for the light cuive of SN I. 

In the present series of papers (I - III) , these numerical results 
are described in detail and the various types of white dwarf models for 
SN I are presented. The preliminary results have been published in the 
proceedings of several workshop and symposia (Nomoto 1980a, b,c; see also 
Woosxey, Weaver, and Taam 1980). In this Paper I, I discuss the presupor- 
nova evolution of the white dwarfs, namely, growth of a helium zone, 
ignition and progress of the strong flash, and the transition from 
hydrostatic through hydrodynamical stages. Cases with the relatively 


•low mv:cr«tion and high initial maaa of tha vdiita dwarf ara invaatigatad 
aapacially in datail, for which tha rasulta by FS and TM wara only pra- 
liminary. Than tha refinad ralationahip batwaan tha triggaring machaniama 
of SN I and tha conditions in tha binary aystam (i.a., tha accration rata 
and the initial maaa of tha whita dwarf) is prasanted. In tha aubaaquant 
papers (II, III), tha hydrodynamic bahavior of tha supamova explosions, 
i.e., the formation and propagation of the detonation wave and deflagra- 
tion wave will be discussed. These models will be compared with the 
observations of SN I. 

In the next section, the assumptions and numerical methods ara 
described. Evolution froc* the onset of accretion through tha early 
hydrodynamical stages are discussed in illl and the dependence of the 
mode of ignition on the parameters are summarized in §IV. In §V, the 
various modes of building up the helium zone, i.e., a rapid accretion of 
helium and the various types of hydrogen shell burning in the accreting 
white dwarfs are discussed. 
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II. PHYSICAL INPUT 


Th« method of cooputatlon and th« physical input in the present 
Investlgstlon are mostly the same an those used by Nc»oto and Sugimoto 
(1977) and thus are briefly summarised here. Hereafter, the normal 
notations (see Appendix in Sugimoto and Nomoto 1980) are used unless 
otherwise stated. 

The evolution from the beginning of the accretion all the way through 
explosion was computed by using a Henyey-type hybrid (in^licit-expliclt) 
method to include both thermal and hydrodynamical equations (see Sugimoto 
1970; Nomoto and Sugimoto 1977; Nariai, Nomoto, and Sugimoto 1900) . 
spherical symmetry was assumed. 

For the shock propagation, the artificial viscosity (Richtmeyer and 
Morton 1967) , Q, defined by 


Q 


1 2 ^_ 2 , 3(v/r) . 2 




3r 


if 


9(v/r) 

3r 


< 0 , 


otherwise , 


( 1 ) 


(e.g. Appenzeller 1970) was used, where v denotes the velocity and the 
pareuneter I was assumed to be 2Ar (Ar denotes the distemce between the 
adjacent zone boundaries) . For convection, the time-dependent mixing 
length theory as formulated by Unno (1967) was applied (see Nomoto and 
Sugimoto 1977 for the expression) . The radiative and conductive opacities 
were calculated by using Iben's (1975) analytic approximations except for 
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th« Cgmpton scjittcrlng op«‘:;ity to include th« offeo*: of «l«ctron» 
d*g«nnr«cy (Saibpson 1961) . 

CogapoalticTn of th« C-fO whit* dwarf was aaaunad to b* ■> 0.5 

whare and Xq danot* th* alMindancaa of carbon and oxygon, raapactlvoly. 
Th* equation of stat* for partially relativistic and partially daganarate 
electrons included electron -positron pairs (Suginoto and N«Boto 1075b) . 

The effects of Coulomb interactions were ta)cen into account by applying 
Van Horn's (1969, 1971) approximation for the liquid state and Salpater's 
(1961) formula for the solid state. The Debye function for specific 
heats was included. The effects of crystal izat ion (Polloc)c and Hansen 
197i) were taken into ac \.»mt. 

For the nuclear reactions, the rates given by Fowler, Caughlan, and 
Zimmerman (1975) , the weak and intermediate screening factors by Graboske 
et al. (1973) and the strong screening factors by Itoh at al. (1977, 1979) 
complemented by Alastuey and Jancovici (1976) were used. For the pycno- 
nuclear regime of carbon burning, the reaction rate by Salpeter and Van 
Horn (1969) under the static lattice approximation was used. The nuclear 
statistical equilibrium composition including electron capture rates was 
calculated ly applying the same data used by Yokoi, Neo, and Nomoto (1979). 

It should be noted here that, for some white dwsurfs with relatively 
slow accretion rates, the interior temperature is much lower than in 
normal stars (see §1X1) so that the applicability of the 3a reaction rate 
by Fowler et al. (1973) must l»e examined. For relatively high tempera- 
tures (T > 8 X Io’k) , the 3a process is a resonant reaction. For lower 
tfcu.^eratures, however, the 3a process can no longer be regarded as a 
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resonant raactlon sc nat a nonraaonant reaction rata should ba uaad 
<CasK'.ron 1959) . Since Fowler at al. (1975) give only the resonant 3cx 
reaction rata, I have calculated the nonraaonant reaction rates for the 
Qt-fa and 'Be-K]i reactions according to the recipe of Clayton (1968) by 
using the nuclear data of Barnes (1980) . The ratios thus obtained between 
the nonraaonant (NR) and the resonant (R) reaction rates are given by 


*^Be+a 


<av> 


NR 


«,v>^ 


Be+a 


( 2 ) 


- 2.84 X exp (33.3 - 50.8 T^‘/’)(Tg'/’ - 0.507)"^, 


and 


r 

a •♦•a 


<ov> 


NR 


<ov>^ 


a+a 


(3) 


- 7.20 X exp (10.7 - 29.1 - 0.907)"^, 


for ®Be+a and a*a teactions, respectively, where Tg ^ T/10®K. Since the 
NR reaction rate has a much weaker temperature dependence than the R rate, 
these ratios exceed unity for very lot> temperatures; in fact, > 1 

for To ^ 0.74 ^u^d r >1 for To ^ 0.28. Therefore, the 3a reaction 
rate by Fowler et al. (.1975) has been multiplied by for 0.28 Te 

^ 0.74 and by r„ *r for To 0.28. Since these NR reaction rates 
were not included by either FS or TM, they did not find the ignition of 
helium at temperatures as low as Tg < 0.5. However, the NR reaction rate 
is greatly enhanced by strong screening at high density even for Tg < 0.5 
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•o that th« avolution of a low tanqparaturo whita dwarf should bs groatly 
dlfforcnt from that found in tha praviouu studlaa (aaa fZIZ ) . 


12 


t 


Ill . RESULTS 


a) Initial Models and Parameters 

The pareuneters £or the initial models of the C-fO white dwarfs are 
summarized in Table 1. Hereafter the superscript (0) denotes the initial 
value just before the onset of accretion, and the initial mass is denoted 
by The 1.08 white dwarf model was constructed from the 

1.08 M« C+C cere in the 7 M© red giant atax of Sugimoto and Manoto (1975a> . 
The core was allowed to cool down for 10®yr after the stripping off of 
the hydrogen-rich envelope. The initial models with larger masses were 
constructed from the 1.08 M« model by increasing the mass artificially 
and again permitting them to cool. These white dwarfs have helium enve- 
lopes as thin as i0~ Mq. 

The evolution of a white dwarf during the accretion phase was 
simulated by adopting the same method as Nomoto and Sugimoto (1977), i.e., 
by assuming that tue helium zone grows in mass at a constant accretion 
rate of dM/dt. The effect of the hydrogen shell flashes was neglected; 
the nuclear energy released by the flashes is mostly radiated away during 
the relatively long interflash accretion phases and the gravitational 
energy release by the accreted helium has much larger effects on the deep 
helium layers (Nomoto and Sugimoto 1977) , I computed six cases A-F for 
different accretion rates of helium, dM/dt, and the initial masses, M^+Q' 
as summarized in Table 1. 
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b) l^^o7.utlon During th« Accretion Pheee 


1) Thermal Effects of Accretion 


The accreted matter g^idually formr a helium zone on the C-fO core 

and thus increase the white dwarf mass, M. As a result, the matter in 

2 

tne interior is compressed to higher densities as seen in Figure 1 , where 
the changes in the density distribution for cases A and C are shown. 

Such compression of the interior by accretion releases gravitational 
energy and increases the interior temperature. At the same time, the 
accreting white dwarf loses energy by emitting photons and neutrinos at 
a rate of Lp^j and L^, respectively. 

These ccsnpetitive effects on the interior temperature are shown by 
the energy ccnservation equation as 


M 




3u 


din T 


I’m 


ain T P dt 


dMy 




'^here u denotes the specific internal energy. In the strongly electron- 
degenerate state (i.e., >> 1 where ij; denotes the chemical potential of 

the electron in units of kT and k denotes the Boltzman constant) , the 
coeficient of the Pd (1/p) work term in equation (4) is as small as 


1 - { 


3 In P VI jH P Hi 2\p 


( 5 ) 


* In Figures 1-9, the white d^varf mass is denoted by 
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wher* the Coulomb intexactiona between electrons and ions are neglected 
for B.tinpllclty and and H denote the mean molecular weight per 

electron and per Ion, and atomic mass unit, respectively (see Van Horn 
1971) . It.is implies that most of the gravitational energy release goes 
into increasing the Fermi energy of the electi.‘ons and only a fraction 
contributes to increasing the internal temperature (Van Horn 1971). 
Therefore the effect of the gravitational energy release is negligible 
for the cooling of a constant mass white dwarf. However, the mass accre- 
tion rele£ises a very large gravitational energy so that its thermal 
effect is significant. 

The compression rate due to accretion, Ap ' (din P/dt)^^, is conve- 
niently expressed as 


Ap r (din p/dt)j^ Ap^^^ + 

Ap^”^ Oln p/81n M)g(dln M/dt) , 
Ap^*^^ - - Oln p/3ln q)M(dln M/dt), 


(6) 


where the mass fraction q H Mj,/M for a given Lagrangian shell decreases 

(M) 

as M increases. The Ap term is due to the increase in density at 
fixed q as a result of the increase in M, while the Ap^*^^ term shows the 
compression as the matter moves inward in q-space. In Figure 2, these 
rates Ap, Ap^*^^., and Ap^*^^ are shown as a function of Mj. for the initial 
stage of Case C (M = 1.28 M® and dM/dt =- 7 x 10"’'°M®yr ’) . Since the 
charige in the density distribution as a function of q is quite homologous. 
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(M) 

Oln p/3ln M)q and thus Xp is almost spacially uniform in tha star as 

(q) 

saen in Figure 2. In contrast, Xp is very large near tha surface of 

tha star because of tha steep density gradient there (see Figure 1} . 

(M) 

The energy release due to the Xp term can be estimated from 
equations (4) and (5) as 




(M) ^ 


kT . (K) 


M - 4 X 10 




Xo (M) 

( - - I ".t) 

^ 10 yr 


(7) 


(m) 

where the uniform distribution of T and Xp is assumed. The compressiot. 
{m) {m) 

rate Xp and thus are proportional to dM/dt. Also these values 

are functions of the white dwarf mass, M. In Figure 3, the stelleu: 
radius, R, and the central density, p^, are given as a function of M for 
the models with the interior temperature of 3 x lO^K. Since din P(./dM 
is larger for larger M, especially near the Chandrasekhar limit, the 
compression rate at the center, q - p^./dt, and the corresponding 

Lg are Increasing functions of M. In Figure 3, these functions are 
shown for dM/dt ■ 7 x 10 *°M*yr“*. 

The energy release by the Xp^*^^ term can also be estimated as 


lJ'J' » . P“rdln P f - 3 X XO-% < 7 <«> 


PlH 


(M) 

which is c«nparable to and has a large effect on the temperature 


neeir the surface. 


( 0 ) 


Since the initial luminosities, L . , of the white dwarfs for 

pn 
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cases A-F are in the range of 0.03 L« - 0.06 the conpreseional heating 

(M) la) 

by (L -f L ^' ) are significant as discussed below. 
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ii) Accretion with Moderate Rates (Cases A and B) 

For cases A and B, the ccnnpressional heating rate is appreciably 
higher than the cooling rate of Lpj^. Then the interior temperature rises 
as seen in Figure 4, where the changes in the temperature distribution 
are shown for several stages. Although the accretion begins from the 
same initial model (M^+o “ l-O® M® at t =■ 0 yr) , the thermal history of 
the interior depends sensitively on the accretion rate. 

For Case A, a steep temperature peak appears near the surface as 
seen in the model at t « 1.7 x 10®yr, because the compressional heating 
is more rapid near the surface than in the deep interior (see Figure 2) 
and because its timescale, , is shorter than that of heat transport. 

During later phases, the rise in the temperature becomes slower due to 
the increasing Lp^ and the inward heat transport. The temperature inver- 
sion persists until helium is ignited at t = 2.6 x 10®yr. 

For Case B, the accretion is ten times slower than in Case A and is slow 
enough for the interior to become almost isothermal due to heat diffusion. 
The temperature inversion in the early stage (t = 1.4 x lO^yr) is negligi- 
bly small and disappears at t = 3 x 10 ^yr. Moreover, more heat is lost 
from the star than in Case A while the same amount of helium has been 
accreted. Accordingly, the interior temperature is lower than in Case A. 
Such a dependence on the accretion rate is clearly seen in Figure 5, 
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where the evolutionary paths of the botton of the accreted helium zone 
(hereafter denoted by subscript He and Indicated by the filled circles) 
and the center of the white dwarf (indicated by the open circles) in the 
density - temperature plane are shown for cases A and B. As the mars of 
the accreted helium, increases, and increase. If compared 

at the stage with the S 2 une Pne* **'He lower for the slower accretion. 

When reaches the value given in Table 1, helium is ignited off- 

(1q) 

center at a density of Pjjg ’ , also given in Table 1. The ignition 
point is defined by the time when the nuclear energy generation rate, 

£.^, is equal to the energy loss rate by neutrinos, radiation; 

i.e., Ejj ■ ev + where denotes the radiative flux. The ignition 

was found to take place when the timescale of the temperature rise due 
to nuclear reaction, = c^T/e^, is almost equal to 10®yr. Here Cp 
denotes the specific heat. The corresponding Ignition line is shown in 
p-T plane of Figure 5. For the slower accretion rate, the ignition of 
helium is delayed to a stage with higher Pyg, i.e., larger AMjjg because 
of the lower interior temperature. It should be noted that helium is 
ignited not at the base of the accreted envelope (i.e., Mj. - 1.080 M«) 
but at the outer shell of “ 1.086 for Case A because of the tempe- 

rature inversion in the helium zone. In contrast, ignition occurs at 
the base of the helium envelope for Case B at t ■ 7.2 x 10 ^yr. 

iii) Slow Accretion (Cases C-F) 

For the cases with the accretion rate as slow as dM/dt - (4-7) x 
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10“* ®hJ*yr”* , th« compressional heating rate, is comparable 

to the photon cooling rate, In Figure 6, the evolutional^ paths of 

the bottom of the helium zone and the center are shown for cases C and F 
starting from the same initial model (t ■ 0 yr) but with different dM/dt. 
For both cases, the accretion is slow enough for the interior to bocone 
almost isothermal . 

In the early stages of Case C, the cooling and the heating are 
almost balanced at 0.08 L®, but later when M > 1.33 M®, the tempe- 

rature increases because of the increasing compression rate of Ap ^ (see 
Figure 3) . For Case F, the photon cooling is always danin^mt over the 
gravitational energy release so that the temperature is decreasing during 
the accretion. 

For Case C, T^g is lower than those in cases A and B so that the 
ignition of helium is delayed to the stage with pjjg » 2.4 x 10 g cm” 
and M *= 1.402 M®. This stage is just before the onset of carbon burning 
at the center. For Case F, in contrast, the central density beccxnes high 
enough to ignite carbon by the pycnonuclear reaction before P|jg reaches 
the ignition line of helium. The structure lines at these ignition 
stages (dashed lines) and the ignition line of carbon given by = 10® 

yr are shown in Figure 6. 

Taam (1980b) suggested that accretion slower than 5 x 10~'°M®yr“' 
leads to the ignition of carbon burning at the center prior to the helium 
ignition almost independently of the initial mass of Mq^q. However, he 
did not include the non-resonant regime of the 3a reaction in his compu- 
tation. Since the reson£Uit reaction rate depends much more sensitively 
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on the tenperature than the non-reeonant reaction rate doea, the resonant 
rate drops sharply for T < 5 x lO^K despite the enhancement due to the 
strong screening. Therefore the temperature of the Ignition line would 
be as high as T = 5 x lO^K for p - 5 x in’g cm~’ if only the resonant 
reaction was taken into account. (If the non-resonant reaction had been 
included in Team's (1980b) computation for his models with da/dt - 5 x 
lO'^^M^yr'^ he would have found the helium ignition to occur before 
carbon ignition as can be seen from his Table 3.) 

Because of the weaker temperature dependence of the non-resonant 
reaction rate and the enhancement r'ue to the strong screening, the igni- 
tion line of helium depends strongly on the density around p = (3-5) x 
lO^g cm~’ as seen in Figure 6. This suggests that the mode of ignition 
(helium or carbon) depends strongly on the initial mass of the white 
dwarf for relatively slow accretion rates. 

In order to investigate this dependence, cases D and E were computed 

(see Table 1) . The evolutionary paths in the (Pne» ^He^ ^^c» "^c^ 

plane, and the structure lines at the ignition stages are shown in Figure 
7. For Case E with dM/dt ■ 4 x 10“’ °M*yr"' , the accretion begins from 
^OO which is smaller than for Cass F. It was then found that 

helium is ignited prior to ceurbon, just the opposite result to that for 
Case F. This is because, for the small Mq 4 .q, the mass of the accreted 

helium, ^ue' become large enough for p^g to reach the ignition line 

of helium before the white dwarf mass becones large enough 1.4 M«) to 
ignite carbon. Similar comparison can also be made between cases C and 
0 for which dM/dt again has the same (but larger) value, namely. 
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7 X . For Cas* C which haa the amaller of 1.28 M«, halium 

ignition occura bafora carbon ignition aa deacribed pravioualy, whila for 
Caaa D (Mc+q “ ^ .35 M«) carbon ignition occura firat. 

Thaae caaaa demonatrate how changes in can affect tha nature 

of tha triggaring machaniam for the supernova explosion when dM/dt is 
fixed. It is found that, even for tha very slow accretion, the carbon 
ignition occurs only for the white dwarfs whose is larger than a 

certain critical mass Mc-fO* 

The critical mass, is mainly a function of dM/dt and is obtained 

in the following way. A crucial factor determining the ignition point is 
the temperature at the bottom of the helium zone and the center, i.e., 

Tjje iuid Tj,, which may depend on dM/dt and the initial models of the white 
dwarf (i.e., most cases, the white dwarf is almost 

isothermal, i.e., T^j^/T^ - (0.85 - 0.95) (see Table 1). Taam (1980b) 
found that the evolutionary paths in the - T^, plane for the models 
with the same dM/dt and do not depend on at later stages. In 

the present invest: jation, models starting from different were ccan- 

puted for the same dM/dt (Table 1). For both cases E and F with dM/dt 
■ 4 X lO'^’M^yr"', the central temperature reaches almost the same value 
of 1.7 X lO^K at later stages with p^. = (0.8-1) x lo'°g cm”^. Also for 
dM/dt - 7 X 10~'®Moyr”', T^, 3.7 x lO^K is reached for both cases C amd 

D when the flash is ignited. 

These results show that T^ emd the interior temperature distribution 
near the stage of carbon ignition depend mainly on dM/dt and only slightly 
on the initial conditions (i.e., M^+q This is because the 
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whit* dwarf avolv** gradually toward* th* *tat* with 

(•** also FS) . Ne iiote that, in th* pr***nt computation, th* t*np*rature 
x£j somewhat low*r than in FS and TM for th* cases with the same dfVdt 
because th* effect of Debye cooling was not included in FS and TM. 

Once the p~T structure line at the stage of carboii ignition is given 
as a function of dM/dt, can be estimated as follows. For dM/dt ■ 

4 X lO'^^M^yr'^ the structure line crosses the ignition line of helium 
at “ 4.5 X 10 ®g cm"^vd»ich is th* density at the shell with ■ 

1.17 M« in th* model just at carbon ignition (Case F) . This implies that 
Mc*o ~ 1>17 M« for this accretion rate. Similarly M^q is estimated to 
be 1.28 for dM/dt ■ 7 x 10' * “M*yr"* . For other values of dM/dt, the 
temperature and the density distribution in the model just at the carbon 
ignition are obtained by interpolation and extrapolation} M^*q is then 
obtained as a function of dM/dt. The line of - M^,*q in Figure 8 

discriminates the mode of ignition in the (M^^q - dM/dt) plane. 

If the pycnonucleer reaction rate of carbon burning under the relaxed 
lattice approximation (Salpeter and Van Horn 1969) is adopted, the 
ignition density of carbon is lower than that based on the static lattice 
approximation by a factor of 1.3 (Figure 6). Correspondingly, the criti- 
cal mass, Mq*q, is smaller as indicated by the dashed line in Figure 8. 

c) Initiation of Thermonuclear Explosion 

As described in the preceeding subsection, accretion onto C'fO white 
dwarfs ignites various types of unstable nuclear bui.ning depending on 
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dM/dt and In this siibaection, I describe the transition phase 

from the ignition through the explosive stages. 

i) Strong Helium Shell Flash 

The progress of the unstable helium shell burning (flash) has been 
computed for cases A-C. Since the hydrostatic stages of the shell flashes 
under similar conditions have been investigated by FS and TM, I will 
focus on the transition from the hydrostatic through hydrodynamic phases 
of the flash. 

Because of strong electron degeneracy, the temperature of the burning 
shell, Tj|g, rises as shown by the solid lines in Figures 5 and 6. A 
convective zone develops, and the superadiabatic temperature gradient 
increases as Cp increases. When Tj^g reaches the stage with Tp ' lOOx^ 

(Tjj » (24TTGp)"'/^ is the dynamical timescale), connection can transport 
only a negligible fraction of the nuclear energy generation. The dashed 
lines in Figure 5 snows the structure lines of the white dwarfs in p-T 
plane when such a blocking of heat occurs. Afterwards the progress of 
the flash is confined to a very thin shell euid rises quickly. 

Soon Tj^g reaches a deflagration temperature , T^^^, defined by Tp ■ 
as indicated in Figures 5 2 md 6. The rapid nuclear energy release 
increases the pressure; the overpressure AP, i.e., the excess of pressure 
over that required for the hydrostatic equilibrium, develops. In other 
words, the hydrodyn£unic phase begins. The effect of the inertial force 
delays the expansion of the burning shell emd keeps T|{g rising (Sugimoto 
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1964), Th« ratio AF/P increases fron 8 x lO'** to 0.1 (Case A), from 4 x 
10"’ to 3 X 10"* (Case B) , and from 5 x 10"® to 7 x 10"* (Casa C) during 
the evolution fr(»> the stage with T„g - T^ef to - 2 the thermal 

overpressure is larger for the lower density because of lower electron- 
degeneracy. Ultimately, the increaaing Ap/p results in the formation of 
a shock front and an off-center detonation wave as will be described in 
Paper II (see preliminary reports by Nomoto (1980a, b)). 

As has been shown by Nomoto and Sugimoto (1977) for the off-center 
flash in helium white dwarfs and by FS and TM for C+0 white dwarfs, the 
accretion at the higher rate considered here (i.e., dM/dt ^ 4 x 10"*M®yr"') 
ignites the helium at such a low density, 10®g cm"*, that the 

resultant helium shell flash is too weak to induce any dynamical effects. 
Therefore, accretion can trigger an off-center helium detonation only in 
a restricted region of the (dM/dt, plane, as shown in Figure 8. 

ii) Carbon Flash 

The progress of the carbon burning for Case F is shown in Figure 6. 

Since the pycnonuclear reaction rate depends very weakly on the tempera- 
ture, the central temperature rises relatively slowly with a timescale 
of dt/dln Tj, = 10®yr. During such a slow rise in Tg, the central density 
continues to increase due to the accretion, which in turn increases the 
reaction rate. For T > '.i x 10*K, the carbon burning changes from the 
pycnonuclear regime to the thermonuclear regime so that the thermal run- 
away is accelerated (see also Taam 1980b) , leading to the carbon deflagration 
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HUparnova (Nomoto «t al. 1976} Nomoto 1980a, b). Tliia type of auparnova 
axplosion occura for Mq^q - in the (dM/dt, Mq 4 .q) plana of Figure 8. 

The carbon deflagration aupemova ie also triggered by rapid accre- 
tion in the region of dM/dt 4 x 10"*M«yr“* in Figure 8. As mentioned 
above, the helium shell flashes for such an accretion are very weak and 
so will recur many times to increase the mass of the C-*-0 white dwarfs 
Such ar> evolution is very similar to the growth of the degenerate C-t-O 
core in red giant stars, although further investigation is needed on the 
mixing of hydrogen-rich material into the helium-burning convective rone 
which might lead to the hydrogen flash and s-process nucleosynthesis and 
make the evolution more complex (FS; TM; Iben 1981) . Finally carbon will 
be ignited at the center of the white dwarf, which will lead to the carbon 
deflagration supernova (Nomoto et al. 1976). 
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IV. CXMCXAJSZON 


B«a«d on th« raaults of tho prscadlng section, the relationship iMtween 
t!te triggering mechanisms of SN I and the conditions in the binary system 
is oummarired in the » dM/dt) plane (Figure 6) . SN I are triggered 

by either the oJf-center detonation of helium for slow and intermediate 
accretion rates or the centrr.l carbon flash for slow and rapid accretion 
rates. Such a variation in the mode of ignition may cause some variations 
among SN I {Barbon et al. 1973; Pskovski 1977; Branch 1981). 

This diagram is different from Team's (1980b) preliminary one in the 
part for slow accretion (dM/dt < 1 x lO'^M^yr”') whore both off-certer 
detonation and carbon deflagration are possible; the domains in the (Mq^O 
- dM/dt) plane corresponding to the two triggering mechanisms are separated 
by the line for M^^^q ■ ^he carbon is ignited by the slow accre- 

tion for the case with M^+q ^ ^C+O' the white dwarf has a thick helium zone 
of mass AMjig - 1.40 M® - M^^+q - 1*40 M® - M^*q. The upper bounds on the 
mass of the helium zone are 0.25 M® and 0.27 Me for the pycnonuclear reac- 
tion of carbon under the static .d relaxed lattice approximation, respec- 
tively. This is in contrast to the case of carbon ignition under rapid 
accretion in which AM(|® should be smaller than 6 x 10~‘*M® in order to avoid 
the occurence of helium detonation (FS) . The presence of such a thick 
helium zone may also 'ause the variations among SN I. The present finding 
implies that the initial mass function of the white dwars in close binary 
systems (Webbink 1980) is important for estimating the frequency of each 
type of supernova. 
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V. DISCUSSIC*! 


In Figure 8,. dM/dt denotes the accretion rate of helium which is a 
result of hydrogen shell burning or a direct transfer of helium from the 
companion star. When the helium is produced by the recurring hydrogen 
shell flashes, M/At is not necessarily the same as the accretion rate oi 
hydrogen-rich matter, dM^j/dt, from the companion star. In the following, 

I will discuss several m^^des of building up of the helium zone. 

a) Direct Transfer of Helium 

During the evolution in some binary systems, it is possible that 
helium is directly transfored from the companion star to the white dwarf : 

If the companion star is a helium star in the mass range of 1.0 - 2.5 
its helium envelope expands to a red giant size when the degenerate C+0 
core (PaczyAski 1971) or O+Ne+Mg core (Nomoto 1900a, b) is formed. Then 
the helium envelope overflows the Roche lobe and accretes onto the white 
dwarf rather rapidly. Also the emission of gravitational waves in the 
compact binary system could drive the transfer of helium from the degene- 
rate helium dwarf companion to the white dwarf. Although the observed 
examples of such compact binary systems as HZ-29 (Faulkner, Flannery, and 
Warner 1972) and G61-29 (Nather, Robinson, and stover 1981) have too small 
a mass for the system to evolve into a SN I , more massive systems could 
exist . 

In such a direct accretion of helium, very rapid accretion is possible. 
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If the accretion rate of helium, dM/dt, is higher thw a certain critical 
rate, (dM/dt) given in equation (9) below, the accreted helium will 
be piled up to form a red giant like envelope after a very weak flash as 
also shown for the rapid accretion of hydrogen-rich matter (Nomoto, Nariai, 
^u^d Sugimoto 1979) . This critical rate corresponds to the growth rate 
of the mass of the C-t-0 core through helium shell burning at the bottom 
of a helium envelope of red-giant size. Applying Uus's (1970) result for 
the helium star of 1.5 - 3.0 M®, this rate is given by 

(dM/dt) = 7.2 X 10"® (M(,^q/M* - 0.60) M«yr"‘ , (9) 

for 0.75 - - 1.38 emd shown in Figure 8. The upper limit for the 

accretion rate of helium is determined by the Eddington limit, (dM/dt) gjj^, 
and IS also shown in Figure 8. 

Formation of an extended helium envelope will lead to the phase of 
double degenerate cores embedded in the car.uiion helium envelope. More 
investigation is needed to clarify whether such a system evolves into 
twin white dwarfs by blowing off the envelope (Taam, Bodenheimer, and 
Ostriker 1978) or into a supernova triggered by the collision of the cores 
(Sparks and Stecher 1974) . If the latter case were to occur, a peculiar 
SN I might result. The outcome for the accretion of dM/dt < (dM/dt) 
have already been summarized in Figure 8. 
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b) Hydrogen Shell Burning 


When the hydrogen -rich matter accretes onto the white dwarf, hydrogen 
shell burning is ignited. How much mass of ac«.. reted hydrogen is processed 
into helium depends on the type and the strength of the hydrogen shell 
burning and thus depends sensitively on the accretion rate of hydrogen- 
rich matter, dM^j/dt. This is analogous to the case of helium shell burning 
discussed in the preceding subsection and summarized as follows (see 
Sugimoto and Miyaji 1980 for a review) . 

(1) Rapid Accretion Forming a Red Gieint Like Envelope ; When the 
accretion is as rapid as dMj^/dt % (dM/dt)jy^, the accreted matter forms a 
red giant like envelope (Nomoto et al . 1979). This critical rate, 

(dM/dt)j^, corresponds to the growth rate of the degenerate core in red 
giant stars due to the hydrogen shell burning. Paczyiliski ' s (1970) core 
mass to luminosity relation assuming a hydrogen abtindance of X =• 0.7 gives 

(dM/dt)jy^ - 8.5 X 10"’ (M/M® - 0.52) M«yr"*, (10) 

for 0.60 5 M/M® 5 1.39, where the core mass is replaced by the white dwarf 
mass M. In Figure 9, (dM/dt)jy^ and the Eddington imit, (dM/dt)gj^, for 
hydrogen are shown. The red giant like envelope f the white dwarf can 
cause various further evolution of the interacting binary as mentioned 
for the rapid accretion of helium; if the common )iydrogen-rich envelope 
is blown off by spiraling of the white dwarf into the ccxnpanion star 
(Taam et al. 1978), a helium companion star would be left, and the resultant 
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system, having a relatively short period, could evolve into a phase of 
direct accretion of helium onto the vdilte dwarf. 

(2) Steady Burning ; For the accretion rate in the range of 0.4 
(dM/dt)|^ % <U4||/dt < (dM/dt)|^, the hydrogen shell burning is stable 
(Paczyi^ski and Zytkow 1976; Sion, Acierno, and TomcciZyk 1979; sienkiewics 
1980) . Such hydrogen burning processes the accreted matter into helium 
at a rate of dM/dt ■ dM{j/dt. 

(3) Recurrence of Flashes ; When the accretion rate is slower than 
0.4 (dM/dt) the hydrogen shell burning becomes unstable emd flashes. 

How much of the accreted mass is processed into helium without escaping 
frcxn the white dwari depends on the strength of the flash and thus on 
dM||/dt. For higher dM||/dt, the mass of the accreted matter, AM||, at the 
time of hydrogen ignition is smaller as shown in Figure 9. Here the 
value of AM|| is taken from the models computed by Narial and Nomoto (1979) 
where the compressional heating due to accretion is just balemced with 
the cooling due to heat conduction; such a condition is realized after 
many flash cycles. 

In general, smaller AMj^ emd smaller M lead to a weaker flash because 
of the lower pressure at the flashing shell (Sugimoto et al. 1979). When 
AM^ is sufficiently small, the flash gives rise to no mass ejection nor 
appreciable expansion so that the white dwarf appears as em EUV nova 
(Shara, Prialnik send Shaviv 1978) . Sion et al. (1979) also found that 
such weak flashes occur in the models with M <■ 1.2 and dM|j/dt >■ 2 x 10~^ 
- 10“®M*yr*'. Therefore, a large portion of the accreted matter will be 
processed into helium for accretion as rapid as dM^/dt = 2 x 10~^ - 
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10“®M*yr”' (see also Nomoto 1980c) . 

For the slower accretion, the flash is stronger because of larger 
AM}{. It is still li)cely that a finite amount of helium is left even when 
a nova-lik‘3 explosion occurs, so that the helium zone can grow at a rate 
of <3M/dt which is slower than dMy/dt. One might worry that the companion 
star could not supply enough hydrogen to lead to SN I if dM/dt is much 
smaller than dM^/dt by, say, a factor of 10. It depends on the parameters 
in the binary system. If the initial white dwarf is as massive as Mq+q 
~ 1.3 - 1.4 Mg, the accumulation of at least 10” ‘ - 10”^Mg helium can 
trigger the helium detonation (PS) . Therefore, the accretion of in total 
0.1 - 0.01 M« hydrogen is enough to give rise to SN I. If the companion 
star would lose all its hydrogen-rich envelope before a SN I is triggered, 
the accretion would stop and the white dwarf would cool again; as a result 
of further evolution, however, it could occur that the ccxnpanion star 
would transfer helium onto the white dwarf driven by the expansion of 
the helium envelope or by the emission of the gravitational waves. 

(4) Diffusion for Slow Accretion ; When the accretion rate is as 
slow as 10”*° - 10”’*Mgyr“*, the accretion continues as long as ' 

Al^/ (dMj^/dt) - 10° - lO^yr before the ignition of the hydrogen shell 
burning as estimated from Figure 9 . Starrfield, Truran, and Sparks (1981) 
suggested that the CNO nuclei could diffuse out of the accreted matter 
in such a long tacc* leads to stable bu’-ning by the proton-proton 

chain, a helium zone will be built up at a rate of dM/dt = dM^/dt. 

However, more quantitative investigation of the diffusion process during 
the accretion is needed. 
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The above discueslon In points (l)-(4) shows that SN X from accreting 
white dwarfs can occur over a rather wide range of accretion rates. If 
other conditions are the same, more rapid accretion triggers SN I more 
frequently because of the shorter accretion time. In this respect, the 
investigation of the connection between SN I and the symbiotic stars might 
be interesting, because Paczynskl and Rudak (1980) suggested that some 
symbiotic stars may consist of white dwarfs on which hydrogen is burning 
steadily or weak flashes are occurring. Such weak hydrogen flashes recur 
many times before a SN I is triggered; the number of flashes is estimated 
from the mass ratio between at the initiation of flash (helium or 

carbon) and the average value of AMj^ during the growth of M assuming 
negligible mass ejection. The value of AM^ gives an upper limit for the 
hydrogen mass when a SN I occurs. Their smallness compared with M is 
consistent with the hydrogen -deficiency in SN I. 

It is also import^ult to Investigate various evolutionary paths of 
.-acting binaries leading to a phase of direct transfer of helium from 
a companion star to the white dwarf (see reviews by Webbink 1979 and by 
Tutukov 1980), e.g., the evolution after the contact system is formed and 
after the recurrent nova explosions finally cease as discussed above in 
point (1) and (3) . 
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TABLE 1 


Initial and Ignition Modala ^ 


Initial Modala 


Case 

A 

B 

C 

D 

E 

F 


i.oe 

1.08 

1.28 

1.35 

1.13 

1.28 

dM/dt lM*yr"^) 

3x10"* 

3x10"* 

7x10"^® 

7xl0‘ * ® 

4x10" ‘® 

4x10"* ® 

(t.) 

0.03 

0.03 

0.07 

0.08 

0.03 

0.07 

log T^0> (K) 

7.37 

7.37 

7.42 

7.42 

7.34 

7.42 

log pjo) (g cro“*) 

7.76 

7.76 

8.59 

9.14 

7.92 

8.59 


Ignition Models 


Case 

A 

B 

C 

D 

E 

F 

Ignited Fuel 

He 

He 

He 

C 

He 

C 

t (yr) 

2.60x10® 

7.77xl0’ 

1.74x10® 

7.57xl0’ 

6.68x10* 

3.03x10 

'^He 

0.078 

0.233 

0.122 

0.053 

0.267 

0.121 

log Tne (K) 

7.90^ 

7.77 

7.53 

7.49 

7.20 

7.16 

log pHe<9 cm"*) 

6.44^^ 

7.56 

8.38 

7.93 

8.67 

8.36 

log Tc (K) 

7.53 

7.80 

7.58 

7,56 

7.23 

7.22 

log Pg(g cm"*) 

8.02 

8.82 

10.04 

10.08 

9.90 

10.03 


a Subscripts He and c denote the bottom of the helium envelope amd the 
center, respectively. 

b For Case A, P|{g and correspond to the values at the ignited shell 
which does not coincide with the bottom of the helium envelope. 
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FIGURE CAPTIONS 


Figure l> Ctumgea in the density distribution in the accreting white 

dwarfs are shown for cases A and C. The interior matter is compressed 
to higher densities as the white dwarf mass, increases. The 

accretion starts at an initial mass of ■ 1.060 (Case A) and 
1.280 (Case C) . Helium is ignited at ■ 1.156 (Case A) 
and 1.402 (Case C) . Open and filled circles indicate the center 
of the white dwarf and thr bottom of the accreted helium zone, 
respectively. 


Figure 2; The rate.*; of compression due to mass accretion onto the white 
dwarf are shown for the initial stage of Case C (M^d ■■ 1.28 M« and 
dM/dt ■ 7 X 10“* ®M*yr"* ) . These rates are difined by = 

Oln p/9ln M) (din M/dt) , = -Oln p/3ln q)j„j (din M/dt) , and 


(M) ^ (,) 


Figure 3: The central density, p^, and the stellar radius, R, for white 

dwarfs with a central temperature Tq ■ 3 x lO^K are shown as a 
function of the white dwarf mass The compression rate by 

accretion, Xp^j, ^ (<jln p^/dM) (dM/dt) and the corresponding gravita- 
tional energy release for dM/dt ■ 7 x lO'^^M^yr”* and T^ 

3 x lO^K are given as a function of M|(p. 
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Plgur* 4: Changas in tha tamparature diatribution for casaa P, and B ara 

shown. Tha Initial modal (t ■ 0 yr; - 1.08 M^) at tha onsat of 

accretion are the same for both cases A and B. At the stages of 
t <■ 2.6 X 10 ^yr for Case A and t « 7.2 x 10 ^yr for Case B, helium 
ignition occurs. 

Figure 5; Evolutionary paths (solid lines) of the density eutd tempara- 
ture at the bottom of the accreted helium zone (filled circle) and 
at ihe center of the vAiite dwarf (open circle) for cases A (dM/dt ■ 

3 X 10"®M*yr"‘) and B (3 x 10'®M«yr“‘) . For Case A, the rise in 
temperature after the helium ignit.'.on is shown for the shell with 
maximum temperature which does not coincide with the bottom of the 
helium layer (see text) . Dashed lines show the structure of the 
white dwarfs from the center to the surface at the stages when the 
accretion begins (t » 0 yr) / nd when the heat blockage occurs in 
the helium-burning shell. Dotted lines are the ignition lines for 
helium and carbon defined by ■ 10®yr and = £v Th® deflagra- 
tion temperature for helium defined by Tjjg = is also shown by 
dotted line. 

Figure 6; Same as Figure 5 but for cases C (dM/dt = 7 x lO'^^M^yr”') 
and F (4 X 10"^ °M«yr”* ) . The structure lines (dashed lines) are 
shown for the initial stage (t * 0 yr) and for the stage when helium 
(Case C) or carbon (Case F) aure ignited. The ignition lines indi- 
cated by (S) euid (R) for carbon burning (dash-dotted line) correspond 


- 41 



to the pycnonuclear reaction under the static and relaxed lattice 
approximation, respectively. 

Figure 7: Same as Figures 4 and 5 but for Cases D and E. 

Figure 8: Relationship between the triggering mecheuiisms for Type I 

supernovae and the conditions in the binary system, i.e., the initial 
mass of the white dwarf, Mq^.q, and the accretion rate of helium, 
dM/dt. Filled circles show cases A-F computed in the present study. 
Solid and dashed lines .for the critical mass, correspond to the 

pycnonuclear reaction of ceurbon under the static and relaxed lattice 
approximation, respectively. (dM/dt) corresponds to the growth 
rate of the C+0 core in the red giant-size helium star, and (dM/dt) gjjg 
shows the Eddington's critical rate for helium. For (dM/dt) pjjg < 
dM/dt < (dM/dt) gfjg, a red giant-like helium envelope is formed. 

Figure 9; Types and strengths of the hydrogen shell burning are shown as 
a function of accretion rate of hydrogen-rich matter, dM^/dt, and 
the white dwarf mass, M^jjj. (dM/dt) gjj is the Eddington's critical 
rate for hydrogen euid (dM/dt) pjj corresponds to the growth rate of the 
degenerate core in red giant stars. For (dM/dt) gjj 5 dMu/dt < (dM/dt) eh» 
the accretion leads to the formation of a red giemt-like envelope. 

For 0.4 (dM/dt) pu ^ dMp/dt < (dM/dt) pp, the resultant hydrogen shell 
burning is st£d3le and steady. For slower accretion, the hydrogen shell 
burning becomes unstable and flashes. AMp is the mass of the accreted 
hydrogen-rich envelope at the ignition of hydrogen (see text) . 
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